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Abstract: Borexino is an organic liquid scintillator detector located in the underground Gran Sasso National Laboratory 
(Italy). It is devoted mainly to the real time spectroscopy of low energy solar neutrinos via the elastic scattering on 
electrons in the target mass. The data taking campaign started in 2007 and led to key measurements of ^Be and **B 
solar neutrinos as well as antineutrinos from the earth (geo-neutrinos) and from nuclear power reactors. Borexino is 
also a powerful tool for the study of cosmic muons that penetrate the Gran Sasso rock coverage and thereby induced 
signals such as neutrons and radioactive isotopes which are today of critical importance for upcoming dark matter and 
neutrino physics experiments. Having reached 4y of continuous data taking we analyze here the muon signal and its 
possible modulation. The muon flux is measured to be (3.41 ± 0.01) ■ 10~ m~^s~^. A modulation of this signal with 
a yearly period is observed with an amplitude of (1.29±0.07)% and a phase of (179±6) d, corresponding to June 28"*. 
Muon rate fluctuations are compared to fluctuations in the atmospheric temperature on a daily base, exploiting the most 
complete atmospheric data and models available. The distributions are shown to be positively correlated and the effective 
temperature coefficient is measured to be qt = 0.93 ± 0.04. This result is in good agreement with the expectations of 
the kaon-inclusive model at the laboratory site and represents an improvement over previous measurements performed at 
the same depth. 

Keywords: Borexino, Muon, Cosmic, Seasonal, Modulation 



o 

rn 

o 



X 



I. Introduction. The flux of cosmic muons detected deep 
underground shows variations which are in first approxima- 
tion seasonal. The effect is known and studied since many 
decades HI. At Gran Sasso National Laboratory in cen- 
tral Italy the rock coverage is about 3800 m w.e. and the 
expected amplitude of the modulation is ~1.5%. Borex- 
ino is a scintillator detector with an active mass for muon 
detection of 1.33 kt and, being spherical, its acceptance is 
independent of the angle of the incoming muons. It there- 
fore plays a key role in measuring the cosmic muon flux 
and its modulation with reduced systematics. Moreover as 
air temperature data is available from weather forecast web 
services, the correlation with the muon flux can be inves- 
tigated and the effective temperature coefficient can be de- 
termined. Such temperature coefficient, with a larger expo- 
sure, can in future be used to determine indirectly the K/tt 
ratio in the interaction of primary cosmic rays in the atmo- 
sphere, probing a complementary energy region compared 
with existing accelerator experiments. 

II. Borexino Detector. The Borexino detector was de- 
signed to have very low intrinsic background. The cen- 
tral scintillation volume, 278 1 of ultra-pure PC (pseu- 
documene) doped with 1.5 g/1 of the fluor PPO (2,5- 
diphenyloxazole), is contained in a spherical Inner Vessel 
(IV), 8.5 m in diameter, made of 125 |j m thick nylon. It is 
shielded by two buffer layers consisting of PC and a small 



amount of the light quencher DMP (dime thy Iphthalate). 
The surrounding Stainless Steel Sphere (SSS) of 13.7 m di- 
ameter holds 2212 inward-facing 8" photomultiplier tubes 
(PMT) that detect scintillation light from the central region. 
All these components form the Inner Detector (ID) ||2|- 

Though Borexino is located deeply underground, in Hall C 
of the Gran Sasso Laboratory (LNGS), the residual muon 
flux is ^1.2 |j/m^/h, still too large for neutrino measure- 
ment, so the muons must be individually tagged. To ac- 
complish this task the ID is suiTounded by a powerful muon 
detector Q . It is composed by a high domed steel tank of 
18 m diameter and 16.9mheight filled with2 lOOtof ultra- 
pure water and instrumented with 208 PMTs which detect 
the muon Cerenkov emission. The Water Tank also serves 
as additional passive shielding against external radiation. 
This system is called the Outer Detector (OD). 

III. Cosmic Muon Flux. This analysis is based on the 
first 4 years of Borexino data, taken between May 16*'' 
2007 and May 15*'' 201 1, with the exclusion of calibration 
data, data not passing the validation procedure and data for 
which the OD was not functioning properly. Events in co- 
incidence with the spills of the Cern-to-GranSasso neutrino 
beam are discarded (details in t3J). The remaining data set 
shows no prolonged or unevenly distributed off time. This 
analysis is based only on muon events that triggered both 
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ID and OD. The total resulting exposure is ^ 1.41 • 10*' t d 
and includes a sample of ^ 4.6 • 10^ muons. 

We have measured the muon rate through the ID using 
different strategies at our disposal and achieved identical 
results. The overall detector's efficiency is 99.992%. In 
|[3l are reported details on the muon tagging methods and 
on how the efficiencies have been evaluated. The aver- 
age muon rate is (4310±10) counts per day, where the 
statistical error is negligible and the systematic error re- 
flects the uncertainty in the efficiency and possible thresh- 
old effects. The rate corresponds to a cosmic muon flux of 
(3.41 ± 0.01) ■ lO^^'m^-^s^^, taking into account also the 
uncertainty in the SSS radius. 

This is the first measurement available for Hall C so far and 
the first obtained with a spherical detector at LNGS: ex- 
isting measurements were obtained with detectors whose 
acceptance strongly depended on the muon incidence an- 
gle and are therefore affected by larger systematics. They 
have been performed by LVD in Hall A (|4j) and by 
MACRO in Hah B ([5]) and are respectively (3.31 ±0.03) • 
10-^m-2s-i and (3.22 ± 0.08) ■ lO-^-^s"!. 

IV. Flux Modulation. Muons observed in underground 
sites arise mostly from the decay of pions and kaons pro- 
duced by primary cosmic ray particles interacting with nu- 
clei in the atmosphere[6). Only mesons decaying before 
further interaction produce muons energetic enough to tra- 
verse the rock coverage of an underground site. Air temper- 
ature increases during summer, leading to an expansion of 
the traversed medium, which in turn increases the fraction 
of such mesons. The following formula is generally used 
to relate the muon intensity variations to the atmospheric 
temperature fluctuations: 



(1) 



where 1° = Ii^i_{To,E > Ethr) is the differential 
muon intensity integrated from the energy threshold {Ethr 
^ 1.8 TeV)!?] to infinity, assuming the atmosphere is 
isothermal at temperature T^, and A/^ are fluctuations 
about a{X) is the temperature coefficient that relates 
fluctuations in the atmospheric temperature at depth X, 
/S.T{X)/T^{X), to the fluctuations in the integral muon 
intensity; the integral extends over atmospheric depth from 
the altitude of muon production to the ground. 

Other underground experiments have studied these effects, 
at the Gran Sasso site (MACRO |8 |, LVD |4l) and at dif- 
ferent underground locations( [9J, |10| andrefs. therein). 

The muon intensity measured day by day is shown in fig. 
[T] (lower panel) for the 1329 days for which valid data was 
available. A modulation is clearly visible; fitting the distri- 
bution with the following function: 
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we obtain an average intensity — (3.414 ± 0.002) ■ 
10~''m^^s^^, a period T = (366 ± 3)d, a modulation am- 



pHtude 61^ = (4.4 ± 0.2) • lO^^m^^s"^ coiTesponding to 
(1.29±0.07)%andaphaseto = (179±6)d; the x^/NDF is 
1558/1325. A Lomb-Scargle analysis of the data identifies 
the same period. It should be noted that due to the lim- 
ited size of the detector, a day with 100% duty cycle fea- 
tures a statistical error of ^^1.5% (la), comparable to the 
expected modulation. Therefore in spite of the fair value 
of reduced x^, we regard this fitting exercise only as a first 
order approximation; in the hypothesis, explored here, that 
the modulation is related to the air temperature fluctuations, 
the main maxima and minima can occur at different dates 
in successive years and short term effects are well expected 
to perturb the overall seasonal behavior 

V. Atmospheric Model. The atmosphere consists of many 
layers that vary continuously in temperature and pressure. 
A possible parametrization ( ifTOl and with more detailsQ) 
considers the atmosphere as an isothermal body with an ef- 
fective temperature, Teff, obtained from a weighted average 
over atmospheric depth: 



j;^ dXT{X)W{X) 
J,°°dXW{X) 



(3) 



where the weight W{X) reflects the temperature depen- 
dence of the production of mesons in the atmosphere and 
their decay into muons that can be observed at depth. 

Eq. [T]can now be written in term of the "effective temper- 
ature coefficient" ax- 
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The weight W{X) can be written as the sum W„ + Wk, 
representing the contribution of pions and kaons to the 
overall variation in muon intensity. 



7 + (7 + l)Bl j,K{X){{Eth cos0)/e^,K)^ 

(5) 



where: 



K{X)^ 



(6) 



The parameters A]^ ^ include the amount of inclusive me- 
son production in the forward fragmentation region, masses 
of mesons and muons, and muon spectral index; the input 
values are = 1 and A]^ = 0.38 • r^/TT, where rj^/T^ 
is the K/tt ratio. The parameters ^ reflect the relative 
atmospheric attenuation of mesons; the threshold energy, 
Eth, is the energy required for a muon to survive to a par- 
ticular depth; the attenuation lengths for the cosmic ray pri- 
maries, pions and kaons are Ajv, A^r and Ak respectively 
with 1/A^ ^ = 1 /An — 1 /A^ The muon spectral index 
is given by 7. The meson critical energy, e^r./f , is the meson 
energy for which decay and interaction have an equal prob- 
ability. The value of {Eth cos 9) used here is the median of 
the distribution. The values for these parameters can be 
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Figure 1 : Upper panel: effective temperature Tgff computed day by day using eq. Q and averaging over the 4 daily 
measurements available. Lower panel: Cosmic Muon Signal. The seasonal modulation is evident. Daily binning. 



found in tab. 1 of ifTOl . with the exception of (£'</, cos 6*) 
which is site dependent and is found by MC simulations. 
At LNGS {Eth COS0) = 1.833 TeV according to |7|. The 
dependency of W{X) on Teff is however moderate. 

Since the temperature is measured at discrete atmospheric 
levels Xn, eq. |3]becomes: 



Weight 
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where W^'^ = W^^iXn). 

Fig. 12] shows the temperature in the atmosphere for the 
LNGS site and the weights used in eq. |7] as functions of 
the pressure levels. As it can be seen the higher layers of 
atmosphere are given a higher weight. Here are produced 
most of the muons which are energetic enough to cross the 
rock coverage of an underground site. Muons produced in 
lower level will be in average less energetic and a larger 
fraction of them lies below threshold. 

VI. Temperature Modulation. The temperature data 
was obtained from the European Center for Medium-range 
Weather Forecasts (ECMWF)|ir| which exploits different 
types of observations (e.g. surface, satellite and upper air 
sounding) at many locations around the planet, and uses a 
global atmospheric model to interpolate to a particular lo- 
cation. In our case, the exact LNGS coordinates have been 
used: 13.578E, 42.454N. Atmospheric temperature is pro- 
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Figure 2: Average temperature (solid red line) and nor- 
malized weight W{X) (black dashed line) as a function 
of pressure levels computed at the LNGS site. 



vided by the model at 37 discrete pressure levels in the [1- 
1000]hPa range (1 hPa = 1.019 g/cm^ ), four times a day 
at 00.00 h, 06.00 h, 12.00 h and 18.00 h. Based on this 
data set, the effective temperature Tgff was calculated using 
eq.|7]four times a dajQ 



1. The analysis in IS) and |2 instead used data from the 
air soundings performed by the Areonautica Militare Italiana 
(AM') L12J ~ 1 30km from the lab. This data set is significantly 
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For the 4y period (Teff) = 220.99K. Fig. [T] (upper panel) 
shows (Teff) for each day obtained averaging the 4 avail- 
able measurement and estimating the error from their vari- 
ance. The fit with a function analogous to eq. |2] returns 
Tjl;* = (221.153 ± 0.007)K, amplitude (2.98 ± 0.01)K 
corresponding to 1.35%, period T = (369.2 ± 0.2)d and 
phase (174.0 ± 0.4)d. However here the error bars are 
much smaller and the x^/NDF is very poor confirming that 
the sinusoidal behavior is only a first order approximation. 
Aside from small scale fluctuations, additional winter max- 
ima can be observed which can be ascribed to the known 
meteorological phenomenon of the Sudden Stratospheric 
Warmings (SSW 1 13 1) and whose effect is sometimes com- 
parable in amplitude with the underlying seasonal modula- 
tion. 
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Figure 3: A/^/ (/^) vs. ATeff/ (Teff). Each point is a day. 

VII. Correlation. Observing fig.[T]the correlation between 
fluctuations in atmospheric temperature and cosmic muon 
flux is evident. To quantify such correlation we plotted 
for every day A/^/(/^} vs AT,{{/{T,ft) in fig. |3] Only 
days with duty cycle > 50% have been included for a to- 
tal of 1 165 days. The correlation coefficient (R-value) be- 
tween these two distributions is 0.60 indicating indeed a 
positive correlation. To determine ar, a linear regression 
was performed accounting for error bars on both axes us- 
ing a numerical minimization method. As a result we ob- 
tain ar = 0.93 ± 0.04 with x^/NDF = 1144/1164. We 
have evaluated the systematic error by varying the assump- 
tions adopted in computing the average flux and tempera- 
ture with respect to the available data set and we have found 
that it is small compared to the statistical error This result 
is consistent and features smaller errors when compared to 
ax — 0.91 ± 0.07, the previous measurement by MACRO 
at Gran Sasso [14). 

The predicted value for ar tends asymptotically to unity 
with increasing depth of the site, as a deeper rock cov- 
erage samples a higher portion of the muon energy spec- 
trum. Fig. |4] shows this behavior along with existing mea- 
surements. The method to compute predicted values as a 
function of site depth is detailed in 1 10|, and for LNGS is 
aT — 0.92 ± 0.02 considering muon production from both 
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Figure 4: Predicted values of ar as a function of detector 
depth and the existing measurements at various depths liTOl . 

pions and kaons. The systematic uncertainty was found by 
modifying the input parameters according to their uncer- 
tainties and recalculating. 

With a longer exposure we foresee to measure ar with bet- 
ter precision and open way to indirect determination of the 
K/tt ratio in the interaction of primary cosmic rays in the 
atmosphere with the method detailed in ifTOl [T] and prob- 
ing a complementary energy region compared with existing 
accelerator experiments. 
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incomplete if compared to ECMWF, although it provided us a 
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